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Abstract 

In [5], the low temperature ratchet current in a multilevel system is considered. In this 
note, we give an explicit expression for it and find its numerical value as the number of 
states goes to infinity. 


1 Introduction 

In this note, we compute the stationary ratchet current in the large system size limit. In [3], the 
authors derive a formula for the occupation of a general multilevel system at low temperature. 
As an application, they consider a continuous time version of Parrondo’s game at low temper¬ 
ature (see [3]) and give an expression for the ratchet current. We consider a multilevel system 
determined by a finite number of states. The set of all states is denoted by K. The ratchet 
is modelised by two rings of N states. In the present section, we recall the definitions and re¬ 
sults from [21 section 3] and in the next section we give an explicit expression for the ratchet 
current using the Tutte matrix tree theorem and find its limit as the number of states goes to 
infinity. The states on the outer ring are denoted by (0, i) and on the inner ring by (1, i), where 
i = 1,..., N. The energies are denoted hy Ei, i = 1,... ,N and are such that Ei < ■ ■ ■ < Ejq. 
The transition rates on the outer ring are given by 

A((f, 0), {i + 1,0)) = A((z -b 1,0), (i, 0)) = 

where /3 is the inverse temperature. On the inner ring, the transition rates are constant and 
equal to one, that is, 

A((^,l),(^ + l,l)) = A((^ + l,l),(^,l)) = l. 

The two rings are connected with transition rates constant equal to one, 

A((f, n), (z, 1 — n)) = 1, where n = 0,1. 

The zero-temperature logarithmic limit denoted by ^(x, y) is given by 

(j){x,y) ■■= lim ilogA(x,i/). 

/3->-oo p 

The zero-temperature logarithmic limit of the escape rates of state x is denoted by r(x) and 
given by 

r(a;) := - lim a log (V A(x, y)) = - max(()(x, y). 


*The author acknowledges support from the Swiss NSF grant 200021_132528/1. 
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The logarithmic-asymptotic transition probability is given by e where 

U{x,y) ■■= -r(a;) - 4’ix,y). 

We have U{x,y) > 0 for all x,y £ K. The smaller U{x,y) is, the larger is the probability of 
transition from state x to state y. Hence, the set of preferred successors of x is defined by 

{y £ K \ U{x,y) = 0}. 

When U{x,y) = 0, the probability of transition from a; to j/ is high. Thus we consider the 
directed graph defined by the vertex set K and edge set {{x,y) \ U{x,y) = 0} where {x,y) 
indicates an oriented edge from x to y. The digraph is represented in Figure [T] below. The 
low temperature asymptotic of the stationary occupation is given in the following theorem from 

0 : 



Figure 1: The directed graph . 


Theorem ([31 Theorem 2.1]). There is e > 0 so that as f3 ^ oo, 

p{x) = + 0(e-^^)) 

with 


0(a;) := nunt/(7i) for U{Tx) ■= ^ U{y,y') and 

{v,v')£T% 

Mx) ■■= n a{y,y') = 

TeM{x) (y,y')(^n 

where the last sum runs over all spanning trees minimizing U{Tx) (i-e. T £ M{x) if Q{x) = 
U{l~x)), o,nd a{x,y) are the reaetivities, which are the sub-exponential part of the transition rates 
X{x,y). 
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Here all the reactivities are constant equal to one, a(a:, y) = 1 for all x,y G K. In the present 
case, for all x G K, there exists an in-spanning tree Tx in K^, so that U{Tx) = 0, and therefore 
0(a;) = 0. Let V be the set of states for which r(a;) = 0, it is given hy V = {(1, 0), {i, l),i = 
1,..., N}. We denote / ~ y if / = g + 0{e~^'^) as /3 —>■ oo. For x GV,we have p{x) ~ \M{x)\/Z, 
where |M(x)| is the number of in-spanning trees in . For x ^ D, the stationary distribution is 
exponentially small since from the theorem it is given by p{x) ~ \M{x)\e^^^^^/Z, with r(x) < 0. 
The stationary ratchet current in the clockwise direction is given by 

j((* + l>0),(*,0))-kj'((* + l>l):(bl)), for i = 


where j{x,y) = \{x,y)p{x) - X{y,x)p{y). 

For i = 1, 

Jfl=j((2,0),(l,0))+j((2,l),(l,l)). 

On the outer ring, we have j{{2, 0), (1,0)) = A((2,0), (1,0))p(2,0) — A((l, 0), (2,0))p(l, 0) with 
A((l, 0), (2,0)) ~ 0, A((2,0), (1,0)) = 

p{2,0) ~ ^ l^(^, 0)1 ^-(g2-gi)/3/2^ 

2 2 

so that j((2,0), (1,0)) ~ |M(2,0)|/Z. 

On the inner ring, we have j((2,1), (1,1)) = A((2,1), (1, l))p(2,1) - A((l, 1), (2, l))p(l, 1) with 

A((2,1),(1,1)) = A((1,1),(2,1)) = 1, 


p(2,l):^ 


|M(2,1)| 




|M(1,1)| 


so that j((2,1), (1,1)) ~ 1)| — |M(1,1)|)/Z. The ratchet current is thus given by 

Jr ~ ^(|M(2,0)| + |M(2,1)| - |M(1,1)|). 

Considering converging arborescences, the Laplacian matrix of a directed graph is defined by 
L = D — A where D is the diagonal out-degree matrix and A = {Aij) is the adjacency matrix 
such that Aij is the number of directed edges from i to j. The rows and columns of L are indexed 
by the vertices of the graph. Here, we index it first by the states on the outer ring then the ones 
on the inner ring, that is (1, 0), (2, 0),..., (iV, 0), (1,1), (2,1),..., {N, 1). The Tutte matrix tree 
theorem (see [T]) relates the number of spanning arborescences converging to x in to the 
cofactors of the Laplacian deiLx^y Let x G K. Then for all y G K, 

\M{x)\ = {-ir+ydetLx,y. 

In particular, for y = x, we have |M(a;)| = detL^,. Therefore we have 

Jr ~ (det L{ 2 ,i) + det L( 2 ,o) ~ det L(i 1 )). 


The Laplacian matrix is given by 


L = 


f.A. 

I B 

[ Id 

: c 


where A is the N x N lower triangular matrix given by 

/I 

-1 1 


A = 


V-1 


-1 1 

0 1 / 
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B is the N X N matrix such that all coefRcients are zero except -B(i,o),(i,i) = the matrix Id 
is the N X N identity matrix and C is the following circulant matrix 


/3 

-1 


C = 


-1 

3 


V-1 




■■■ -1 

-1 3/ 


2 Calculation of the ratchet current 


From [3] , the numerator of Jr is given by 


det F( 2 ,i) + det -^( 2 , 0 ) — det F(i,i) = det Bat-i — 2 det i?Ar _2 — 2 


where is the N x N tridiagonal matrix with 3 on the diagonal and — 1 on the two off-diagonals 
which satisfies the recurrence relation det Bat = 3deti3Ar_i — deti?Ar-2 with detBi = 3 and 
det i ?2 = 8. By solving the associated characteristic equation, it comes 


det Bpf = 


5-3^5^- 

10 V 2 ) 


5-^375/3-hy5\^ 

10 V 2 ) 


The normalisation factor is given by 


xeK Tx (y,z)GTx x&V xeTi 


The sum is over the states in B since the contribution of the states which are not in B is 
exponentially damped. Therefore we have 

N 

Z ~ det Q) + ^ det (1) 

i=l 


We have 


detL(i_o) = detC. 

The circulant matrix C has eigenvalues given by /ij = 3 — 2 cos(27rj/iV), j = 0,1,..., iV — 1 (see 
[2]). Hence 

Af-l 


detL(i.o) = n (3 - 2cos(27rj77V)) = Ul_^{Vb/2) 
j=o 

where Un is the Chebyshev polynomial of the second kind. Thus 


det T(i 0 ) — ^ 


3-^75 


N 


3- V5 


N 


- 2 . 


( 2 ) 


From the Tutte matrix tree theorem, the cofactor (—det L(i_i) is equal to the number of 
converging arborescences to (i, 1) and is equal to the cofactor of the Laplacian where row (i, 1) 
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and any column is removed. Since the only non-zero element of B is in column indexed by (1,1), 
we choose to remove that one, so that 

|M(*,1)| = (-l)('^+*)+('^+i)detL(,,i),(i,i) = (-l)*+MetC(,-i),(i,i) (3) 

since A is lower triangular. On the other hand, by adding to the first column of C all the other 
ones, we have 


det C = 


1 -1 
1 3 

-1 

1 


-1 


•• -1 

-1 3 


N 


— ^(“1)*"'"^ det 


(4) 


i=l 


Putting equations o, ©, m and o together, we have 


Z ~ 2 det C = 2 




N 


3-^5 


N 


- 4. 


Up to exponentially small corrections e the ratchet current is given for all N by 


Jr — 


5 + 3V5 f3 + V5\ 5-3V5 f3-V5\ 5_+_^/3_+j^ 


10 V 2 

5_3^/3_y5^^-2 

5 V 2 


10 V 2 


5 V 2 


N-2 


- 2 j y (^2((3 + V5)/2f + 2((3 - V5)/2)^ - 4^ . 


As a consequence, in the large system size limit the current saturates and has the following limit 


lim Jr ~ - - ■=. 

N-fco 2 1/5 
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